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PACS 03 . 75 . Kk - Dynamic properties of condensates; collective and hydrodynamic excitations, 
superfluid flow 

PACS 03 . 75 . Lm - Tunneling, Josephson effect, BoseEinstein condensates in periodic potentials, 

solitons, vortices, and topological excitations 
PACS 67 . 85 . H j - BoseEinstein condensates in optical potentials 

Abstract. - The phenomenon of dynamical localization of matter wave solitons in optical lattices 
is first demonstrated and the conditions for its existence are discussed. In addition to the trapping 
linear periodic potential we use a periodic modulation of the nonlinearity in space to eliminate 
nonexistence regions of gap-solitons in reciprocal space. We show that when this condition is 
achieved, the observation of dynamical localization in true nonlinear regime becomes possible. The 
results apply to all systems described by the periodic nonlinear Schrodinger equation, including 
Bose- Einstein condensates of ultracold atoms trapped in optical lattices and arrays of waveguides 
or photonic crystals in nonlinear optics. 



An interesting phenomenon occurring in periodic non- 
linear systems is the possibility to sustain different types 
of localized waves for arbitrarily long times due to the in- 
terplay between periodicity, dispersion and nonlinearity. 
In particular, optical lattices (OL) can support so called 
gap-solitons, i.e. localized solutions of the periodic non- 
linear system with energies located inside forbidden gaps 
of the spectrum. These solutions exist both for focus- 
ing (self-attractive) and defocussing (self-repulsing) inter- 
actions due to the coupling induced by the nonlinearity 
between forward and backward propagating waves under- 
going Bragg scattering by the periodic structure [1]. Gap 
solitons can be found in many physical systems includ- 
ing photonic crystals and layered structures [2] , arrays of 
optical fibers [3] and Bose-Einstein condensates (BECs) 
in OL. In this last case, gap-solitons have been experi- 
mentally observed in [4] for BEC with positive scattering 
lengths (repulsive interactions) and their properties have 
been characterized both for discrete (BEC arrays) [5, 6] 
and for continuous BEC [7,8] (see also reviews [9-11]). 
As stationary states with energies inside band gaps they 
are intrinsically localized and remain stable for long time. 



By applying external forces, such as the one generated by 
the gravity (for BEC in periodic vertical traps) or by the 
acceleration of the OL, solitons can be forced to move. A 
moving soliton however has a frequency belonging to an 
allowed band. It is not any more an exact solutions of the 
one-dimensional Gross-Pitaevskii equation with a periodic 
potential, but can be viewed as a persisting wave packet. 

On the other hand, it is known from the linear theory, 
that a time periodic force applied to a quantum particle in 
a lattice can result in the phenomenon of dynamical local- 
ization [12] . In the simplest terms this effect can be viewed 
as a resonance which occurs when the strength of the force 
and its frequency take certain discrete values. Then a 
particle undergoes oscillatory motion in a finite spatial 
domain (namely in this sense the localization is under- 
stood). Outside of these resonances the particle acquires 
a drift velocity which enables transport. For sinusoidal 
time dependent forces resonances are found to be roots 
of the zero-order Bessel function. On a pure quantum 
mechanical level [13] dynamical localization can also be 
characterized as a complete suppression of the inter-well 
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tunneling 1 of the particle in a periodic potential [14, 15]. 

The linear dynamical localization has been experimen- 
tally observed in several systems including arrays of curved 
optical waveguides [16] and BEC trapped [15,17]. Prop- 
erties of the ac forces (electrical fields) that result in the 
dynamical localization of an electron in a periodic poten- 
tial were theoretically investigated in [18] where it was 
shown that for generic (i.e. beyond tight-binding) band 
structures exact dynamical localization is possible only 
for electric field displaying discontinuities at all changes 
of the sign of the field (the dynamical localization phe- 
nomenon, however, was shown to be surprisingly tolerant 
to the smoothing of these discontinuities [19]). 

Thus, while the concept of a gap-soliton is intrinsically 
nonlinear the concept of dynamical localization is inde- 
pendent on nonlinearity. At the same time, the arguments 
applicable to a single particle do not mean of course that 
a wave packet composed of oscillating atoms will oscillate 
without dispersion, i.e. preserving its width. Then the 
natural question arises: Is it possible to observe dynam- 
ical localization of a non spreading wave packets, i.e. of 
gap-solitons? Possibility of positive answer to this ques- 
tion stems from the results reported in Ref [20], where it 
was shown that properly designed nonlinearity can guar- 
antee long-lived Bloch oscillations of gap-solitons of small 
amplitude originated by constant forces. 

The aim of the present Letter is to show for the first 
time how the dynamical localization of gap-solitons of the 
continuous periodic nonlinear Schrodingcr (NLS) models 
can be achieved by means of suitable spatial modulations 
of the nonlinearity (nonlinear OL). For BEC in OL this 
modulation can be obtained by changing the scattering 
length along the condensate using the optically induced 
Feshbach resonance technique. We shall use the spatial 
modulation to suppress (or strongly reduce) non existence 
regions for gap-solitons in the BZ. As a result we have 
that the resonances for the dynamical localization of gap- 
solitons are the same as the one obtained in the seemingly 
very different linear system. 

An intuitive understanding of why the oscillatory mo- 
tion of single particles is accompanied by the spreading 
out of the wave packets and why this situation can be in- 
verted in presence of space dependent nonlinearities can 
be obtained from the analysis of nonlinear discrete models. 
It is known that the usual (e.g. physical) discretization of 
the NLS equation (DNLS) [21] 

ic n = OJ C n + U>i(c n+ i + Cn-i) + x\c n \ 2 C n (1) 

has the onsitc nonlinearity and does not display any dy- 
namical localization phenomenon in presence of ac linear 
forces (any initially localized wave packet quickly spread 
out under the action of the force). In contrast with its con- 
tinuous limit, Eq. (1) is nonintegrable, having only two 
conserved quantities: the energy and the norm. In BEC 

1 The complete suppression of inter- well tunneling is a valid state- 
ment only in the high-frequency limit 
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Fig. 1: (a) The lowest linear band energy E vs crystal momen- 
tum Q (bold line) and the evolution of energy of the soliton 
with norm TV = 5.42 during its movement under the constant 
external force *y(t) = V = 0.001 (thin line). The evolution of 
the soliton crystal momentum is Q(t) — —1 + Ft; (b) The drift 
velocity as a function of the initial crystal momentum of the 
gap-soliton for the nonresonant frequency = 0.001437. The 
linear and nonlinear OLs are fixed as U(x) — — 3cos(2a;) and 
Q(x) = -0.777 + cos(2z). 

contexts this equation arises as a tight binding model of 
the NLS equation and is used to describe arrays of cou- 
pled BECs [5,6]. The integrable discrete version of the 
NLS equation, known as the Ablowitz-Ladik (AL) model, 
has the nonlinearity non local i.e. splitted on two adjacent 
sites: x|c„| 2 (c„_i + c n +i), (which must be substituted in 
(1) instead of |c„| 2 c„) and it is known that it does dis- 
play the dynamical localization phenomenon. The first 
observation of dynamical localization for discrete solitons 
was indeed achieved for this model [22,23]. The general- 
ization of the theory of dynamical localization by includ- 
ing the one-site nonlinearity, as well the understanding 
of the physical picture, was suggested in [23] using the 
discrete nonlinear Schrodinger (DNLS) equation as a gov- 
erning model, where again the same frequencies leading to 
the dynamical localization were obtained. 

To understand why the dynamical localization phe- 
nomenon of undistorted spatially localized wave packets 
occurs in the AL model and does not occur in the DNLS 
chain it is of interest to consider how the stability of lin- 
ear modes of the underlying discrete linear chains (they 
are discrete analogies of the Bloch states of the continuum 
models) is affected by the nonlinearity in the two cases. 
One finds that the stationary solutions of Eq. (1) at differ- 
ent band edges are connected by a staggering transforma- 
tion. It is indeed easy to prove that if c n = exp(— iu>t)f n 
is a stationary solution of Eq.(l) for frequency u> and non- 
linearity x, then c n = exp(— i(2u>o — u>)t — iim)f n is a solu- 
tion for frequency 2wo — ui and nonlinearity x — ~ X- This 
means that if at one edge of the band a mode is mod- 
ulationally unstable (and, hence the systems admits the 
existence of small-amplitude solitons), at the other edge 
of the band for the same nonlinearity the same mode is 
modulationally stable and small amplitude solitons can- 
not exist. For the AL model, however, one finds that to a 
stationary solution c„ for frequency u> and nonlinearity x 
corresponds a solution c„ for frequency 2w — uj and the 
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same nonlinearity X- This means that for the AL model 
it is possible to have stable small-amplitude solitons at 
both edges of the band for a fixed nonlinearity. Actually 
one can prove that in this case solitons can exist for all 
values of the crystal momentum since the corresponding 
plane waves are modulationally unstable in the whole BZ. 
This means that AL solitons can slide along the whole 
band under the action of a time dependent force without 
undergoing any passage through nonexistence region. 

While the AL model itself has restricted physical appli- 
cations, it shows the route for achieving the stability of 
gap-solitons in the whole band with the help of nonlocal 
nonlinear terms. At the same time, such terms naturally 
appear in systems with inhomogeneous nonlincarities (for 
a number of nonlocal Hamiltonian lattices and their prop- 
erties see e.g. [24]). Thus one can expect that the dy- 
namical localization of continuous gap-solitons could be 
achieved in presence of a cc-dependent nonlinearity. This 
conclusion can be reached also from another perspective 
using the effective mass picture. As mentioned before, a 
negative product of the effective mass and effective nonlin- 
earity is required for all points of the BZ to allow the soli- 
ton to move along the whole band (see example in Fig. la). 
While the sign of the effective mass depends on the band 
and therefore on the crystal momentum, the sign of the ef- 
fective nonlinearity can be changed only if the nonlinearity 
is space dependent. For a constant nonlinearity, indeed, 
the soliton will not exist at both edges of the band, this 
leading to its spreading and destruction (this is especially 
true when the non existence regions are sufficiently wide 
with respect to the spectral width of the soliton). 2 

The main idea for dynamical localization in the nonlin- 
ear case then is to suppress (or strongly reduce) non exis- 
tence regions to allow the gap-soliton to safely slide along 
the whole band under the action of the external force. To 
show how this can be done we consider the specific case 
of the NLS equation with periodic potential and ac linear 
force 



iipt = —ip xx + U(x)tp + j(t)xip + Q(x)\ip\ 2 %lj. 



(2) 



In the BEC context the spatial and temporal variables in 
Eq. (2) are measured in the units of d/ir and fr/En, re- 
spectively, while the energy is measured in units of the 
recoil energy En = h 2 Tr 2 /(2md 2 ), where d is the lattice 
period and m is the atomic mass. In a typical situations, 
the external force j(t) arises from the acceleration of the 
linear OL. Then Eq. (2) describes the condensate in the 
moving frame. This implies that the nonlinear OL should 
be also accelerated in the same manner as the linear one. 
Notice, that Eq.(2) appears both in BEC and in nonlinear 
optics contexts although the physical contexts are com- 
pletely different. 

The external force will be assumed to oscillate periodi- 
cally in time. While this is not important for the theory, to 




Fig. 2: Soliton center position X vs time t (upper row, pan- 
els a-c) and particle density \ip\ 2 vs time t and coordinate x 
(middle and lower rows, panels d-i) for linear and nonlinear 
OLs U{x) = -3cos(2z) and Q(x) = -0.777 + cos(2x), and ex- 
ternal force parameters V = —0.001, to — 500.0 (upper and 
middle rows, panels a-f), to = 1000.0 (lower row, panels g- 
i), CI » 0.001306 (left column, panels a,d,g), Q ta 0.001176 
(middle column, panels b,e,h) and fi « 0.001437 (lower panel, 
panels c,f,i). The initial condition at t = is a stationary gap- 
soliton located at X(0) = with energy E = —0.7324 near the 
upper edge of the band (Q(0) = 1.0). In panels a-c bold (red) 
lines represent calculation according to the approximation (4), 
and thin (blue) lines - direct numerical integration of Eq. (2). 



simplify the formulas, and to bring the model closer to the 
experimentally feasible situations in the theory of matter 
waves, we consider the cos-like law 7(f) = T cos [f2(t — to)] 
with amplitude T and frequency fi (time to characterizes 
the initial phase). 

Let us now briefly recall the semiclassical arguments 
which allow one to account for the linear force (see 
e.g. [13]). Identifying the center of a Bloch wave packet, 
X, with the coordinate of a quasiparticle (which in the 
nonlinear case is nothing but a matter wave soliton) one 
finds that it is described by the equations 



X = v(Q) 



dE(q) 



dq 



-7(*), 



(3) 



q=Q 



2 the conditions for optimal design of the nonlinear OL to achieve 
the conditions of the soliton existence in the whole band was depicted 
in Fig. la of Ref. [20]. 



where we also introduce the wave-packet center Q in the 
reciprocal space, use v(Q) for the Bloch velocity and E(Q) 
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for the linear OL energy band structure (see Fig. la). 

Taking into account cos-like dependence of 7(t), the ex- 
pression for the coordinate can be represented as 



X(t) =X - 2W17T 



ttQo - sin[fi(r - to)] j dr, 

(4) 

where we have left only the two leading terms in the 
Fourier expansion of the energy E(q) ss loq + 2u>i cos(7rg), 
and use the notations X(t a ) = X and Q(t ) = Q . In 
this case after one period T = 27r/f2 of the external force 
oscillation, the coordinate of the wavepacket will be 



X[l„ + T) = A-„ + r i; ((;„S./ (1 ( ^ 



T 



(5) 



where Jo{z) is the zero-order Bessel function of the first 
kind and vb = — 2ttu>i sin(7rQo) is the Bloch velocity. 
Thus, the mean velocity of the soliton is given by 



{V s ) = v B (Q ) J 



IT 



(6) 



i.e. it coincides with the Bloch velocity at the time to of 
switching on the time dependent force but scaled by the 
factor Jo (if)- 

As it can be seen from Eq.(6), if O = Q n = Tir/z n 
(z n is the n-th zero of the Bessel function Jo(z n ) = 0), 
then (Vs) — so that the motion becomes perfectly pe- 
riodic X(to + T) = X(to) and the soliton wavepacket is 
dynamically localized. In Fig. lb we depict the soliton 
mean velocity as a function of Q(to) for fixed parameters 
of the OLs and fixed frequency Q. From this figure it is 
evident that the soliton becomes dynamically localized for 
any O at the points Q(to) = to, to = — 1, 0, 1 in the BZ. 
The absence of transport at these points, however, is not 
due to force resonances but simply to the Bloch velocity 
becoming zero (from this point of view the dynamics at 
these points should not be considered as dynamical local- 
ization). 

To check these predictions, we performed direct numer- 
ical integrations of the NLS equation (2) to obtain the 
evolution of the center of the soliton X(t) when an exter- 
nal force starts to oscillate at time t = to 



l(t) = 



Tcos [Q(t 



t < t , 

to)] , t > to, 



(7) 



The results are depicted in Fig. 2. As it can be seen 
from the panel (a), when i7 = Oi (corresponds to the first 
zero of Bessel function), the dependence, calculated from 
approximation (4), exhibits the periodic motion. Never- 
theless, the period of this motion is two times less than 
that of external force oscillations 2n/Cli ks 4810, because 
the period is imposed mainly by the Bloch oscillations, 
whose period is 2/T = 2000.0. The dependencies in up- 
per panel, calculated from direct simulation, exhibit move- 
ment, slightly different from periodic due to the fact, that 
we restricted expansion of E{q) into Fourier series by 
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Fig. 3: Particle density \ip\ 2 vs time t and coordinate x for 
linear OL U(x) — — 3cos(2s), constant nonlinearity Q(x) = 1, 
and external force parameters Y = —0.001, Q ~ 0.001306, to — 
500.0 (panel a), to = 1000.0 (panel b). The initial condition 
at t = is a stationary gap-soliton located at X(Q) = with 
energy E = —0.7324 near the upper edge of the band (Q(0) = 
1.0). 



zeroth and first harmonics only. In cases, when dif- 
fers from fl n (in panel (b) £1 = 0.9f2i and in panel (c) 
fl = l.lfii), the soliton movement looses its periodicity. 

In panels d-i of Fig. 2 we depict the time evolution of 
the matter density as obtained from direct integration of 
the NLS equation. The panels d,e,f correspond to panels 
a,b,c, respectively. Notice that in all the cases the soliton 
is long lived as a consequence of a proper design of the 
OLs. In the nonresonant case Jo (if) (panels e,f) 
the mean velocity of the soliton is (Vs) ~ 0.03722 (panel 
e) and (Vs) « —0.04192, this agreeing very well with the 
value predicted by Eq. (6) (see Fig. 1). 

In panels g-i of Fig. 2 we have shown that localiza- 
tion can also be induced in nonresonant case if the initial 
Bloch velocity at time to is zero, in agreement with what 
predicted by Eq. (6) (notice that to = 1000 and the ex- 
ternal force oscillation starts when the soliton is at the 
bottom of the band Qo = 0). In this case the dynamical 
localization can occur for any frequency Q of the external 
force, as expected from Eq.(5). 

To demonstrate the importance of the nonlinear OL for 
the stability of dynamical localization of solitons, we dis- 
play in Fig. 3a, b, the dynamics of soliton under the pe- 
riodic external force with the same parameters, as those 
of Fig. 2d, e, except for the case of coordinate-independent 
nonlinearity Q(x) = 1. As evident from Fig. 3, in the case 
of constant nonlinearity solitonic oscillations in real space 
are not stable. This happens due the above-described 
mechanism of soliton distortion in the modulationally sta- 
ble regions of BZ. This situation differs considerably from 
the case, depicted in Fig. 2, where special shape of non- 
linear OL results in narrow width of the modulationally 
stable regions in the BZ (see Ref. [20]) so that they do 
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Fig. 4: The same as Fig.2, but for Q « 0.000569 (left column, 
panels a,d,g), Q m 0.000512 (middle column, panels b,e,h) and 
f2 ~ 0.000626 (lower panel, panels c,f,i). 

not affect considerably on the soliton stability during its 
propagation in the real and reciprocal spaces. 

The phenomenon of dynamical localization can occur 
also, when external force oscillation frequency corresponds 
to the higher roots of Bessel function. In other words, 
within the framework of the meanfield theory, there exist 
no lower limit on the frequency leading to the dynami- 
cal localization. Using the well known asymptotic of the 
Bessel functions one obtains the behavior of the low fre- 
quency limit (recall that they are measured in the units 
E r /H): f2 ~ r/ (n — j). Reducing the frequency however, 
generally speaking leads to larger amplitudes of soliton os- 
cillations and more sophisticated trajectories. An example 
of this is reported in Fig. 4 where the dynamical localiza- 
tion is shown for the case = 02- In this case, however, 
we see a small discrepancy between a tight-binding ap- 
proximation and direct numerical simulation of the NLS 
system in Eq.(2) (Fig.4a). 

In closing this letter we wish to remark that periodically 
driven OLs have also been recently used experimentally to 
achieve coherent control of matter waves [25] and theoret- 
ically to investigate phase transition from superfluid to a 
Mott insulator by means of Bose-Hubbard model [14] . In 
both cases an effective renormalization of the tunneling 
matrix element in the Bose-Hubbard hamiltonian was as- 



sumed. In the present Letter we have dealt with a quite 
different situation in which each lattice site is populated 
by many atoms so that the mean field approach is perfectly 
valid. The presence of an inhomogeneous nonlinearity in 
our case permits the existence of gap soliton for all points 
in the BZ. The gap-soliton dynamics is then well described 
by the semiclassical equation of motion and is very stable, 
without any visible increase of the local density, this rul- 
ing out any possibility for occurrence of phase transitions 
or failure of the ID approximation (the latter occurring 
when the two-body interaction energy become compara- 
ble to the transverse kinetic energy) . 

In conclusion, we have shown for the first time the pos- 
sibility of nonlinear dynamical localization of matter wave 
solitons by means of a spatial modulation of the nonlinear- 
ity. In a BEC this can be achieved by changing the scatter- 
ing length by means of optically (or magnetically) induced 
Feshbach resonances while the ac force can be generated 
by a time periodically accelerated OL. For other possi- 
ble applications of the theory, which are arrays of coupled 
nonlinear waveguides, the spatial change of the nonlinear- 
ity can be made by changing the refractive index in the 
fibers along the array with dopants, while the periodic lin- 
ear force could be applied by curving the fibers along the 
propagation direction. 

The phenomenon of dynamical localization of solitons 
in OL is expected then to be observed both in BEC and 
in nonlinear optical systems. 
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